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Abstract 

We consider the TASEP on TL with two blocks of particles having 
different jump rates. We study the large time behavior of particles' 
positions. It depends both on the jump rates and the region we focus 
on, and we determine the complete process diagram. In particular, we 
discover a new transition process in the region where the influence of 
the random and deterministic parts of the initial condition interact. 

Slow particles may create a shock, where the particle density is 
discontinuous and the distribution of a particle's position is asymp- 
totically singular. We determine the diffusion coefficient of the shock 
without using second class particles. 

We also analyze the case where particles are effectively blocked by 
a wall moving with speed equal to their intrinsic jump rate. 

1 Introduction 

We consider the totally asymmetric simple exclusion process (TASEP) on 
Z. This is one of the basic one-dimensional interacting stochastic particle 
systems that, despite its simplicity, is full of interesting features. It consists 
of particles moving to the right by jumps of length one. The jumps happen 
at a given rate (the clocks of different particles are independent), and the 
particles are subject to the exclusion constraint — no site can be occupied by 
more than one particle. This model can also be seen as a growing interface 
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with gradient given by the particle density; it belongs to the KPZ (Kardar- 
Parisi-Zhang) universality class of growth models. Recently the fluctuation 
properties of the TASEP and related models have been studied extensively 
using the techniques from random matrix theory [20,28,31,32]. See also 
[13,36,37] for more recent developments on the case where particles can hop 
in both directions. 

In previous works [4, 6, 8, 30] we analyzed the large time t behavior of 
particles' positions for some non-random initial conditions and uniform jump 
rate (say equal to one). For example, if particles start from 2Z_ then the 
large time macroscopic density is given by 



For large time t, the correlation length scales as t^/^ and the fluctuations 
scale as t^^^. Under an appropriate scaling limit, the joint distributions of 
particles' positions are governed by a process which depends on the density 
gradient: (a) for ^ < 0, it is the Airyi process, (b) for G (0, 1) it is the 
Airy2 process, and at ^ = it is the Airy2^i transition process, see [8]. 
Similarly, one can consider the joint distributions of the current at different 
positions instead of the positions of different particles — the limit processes 
are unchanged. 

In this paper we consider a small variation of the above situation, which 
however shows a number of new nontrivial phenomena. Instead of setting 
the jump rate to 1 for all the particles, we modify the jump rate of the flrst 
M particles and set it equal to a > 0. 

There are a few cases to consider. For example, for < a < 1/2 the 
flrst (slow) particles generate a shock where the macroscopic particle density 
changes discontinuously from 1/2 to 1 — a. The fluctuations on the left of 
the shock are Airyi-distributed on the t^^^ scale, while inside the jam region 
they are GUE(M)-distributed on the t^^^ (diffusion) scale, see the body of 
the paper for detail^. Also, the distribution of a particle's position in the 
shock region has a singularity. 

When a reaches 1/2, the macroscopic density becomes constant (it is 
equal to 1/2 everywhere), but the fluctuations are different. In the simplest 
case of M = 1 slow particle, by Burke's theorem [11], our initial condition is 
equivalent to the deterministic one on Z_ (even sites are occupied as before) 
and to the product of Bernoulli measures with density 1/2 on Z^. Thus, 

^Here GUE(M) stands for the Gaussian Unitary Ensemble oi M x M random matrices. 




(1.1) 
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for a = 1/2 there is a transition region where the influence of the initial 
randomness becomes relevant, but it can not be seen macroscopically. 

On Figure [2] we present the whole process diagram^. One of the goals of 
the present paper is to derive the large time fluctuations' behavior in all of 
its regions. 

Another situation we consider is M = oo and a = 2. Under t^/^-scaling, 
the speed a particles effectively act as a wall moving with speed 1. The 
following n "normal" particles then become like Brownian motions with the 
first one being reflected off the wall and the following ones being reflected off 
each other. The large time fluctuations are then given by the antisymmetric 
GUE(M) process (for fixed time it was characterized in [22], see also [14], 
[15]). This is also closely related to the asymptotics of a certain Markovian 
dynamics for two-dimensional interlacing particle systems with a wall, see 
Section 2.3 of [38] and [9]. Using the relation between last passage directed 
percolation with exp(l) random variables and TASEP, one can predict that 
there should be a relation between the maximum process for the largest 
eigenvalue of the Dyson's Brownian Motion and systems of nonintersecting 
paths with a wall. This relation will be made more precise in [7]. 

Our arguments are based on deriving suitable determinantal expressions 
for the quantities of interest and analyzing the resulting (Fredholm) deter- 
minants asymptotically. In most mathematically rigorous argument 
of that kind would require the evaluation of the asymptotics of the kernel 
under the determinant, as well as some control over the decay of the kernel 
at infinity. This last part is often viewed as a technicality, and we omit tail 
estimates in the present paper. 

In this determinantal approach, the main difficulty typically lies in de- 
riving an integral representation for the kernel before the limit transition; 
evaluating the asymptotics is often quite straightforward via the standard 
steepest descent analysis. However, in the shock case mentioned above, we 
faced a new effect — in the large time limit the kernel diverged. We had been 
puzzled by this difficulty for a while, and we view finding the modification 
of the kernel that solved the problem as our main technical novelty. 

Outline. The rest of the paper is organized as follows. In Section [2] we 
explain the macroscopic picture and describe the process diagram. Then we 
state the results for the different parts of the diagram, which are proven in 
Sections HJIHI In Section [H] we obtain the determinantal correlation struc- 
ture and the associated kernel with a couple of specializations. Finally, in 
Section [7] we consider the reflecting wall situation. 

bit like a phase diagram, but in our case instead of phases and phase transition we 
have Umit processes and transition processes. 
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2 Model and Results 

The continuous time TASEP on Z is a model of interacting particle systems 
in which at every instant at most one particle occupy a site in Z. Particles 
jump by 1 to the right with a given jump rate provided the arriving site is 
empty. As a consequence, particles do not overtake each other. Hence, we 
can assign labels to particles, say particle n has position at time t equal to 
Xn(t). We also denote by f„ the jump rate of particle n. Our convention is to 
consider particles labeled from right to left, i.e., Xnit) > Xn+iit) for any time 
t. We denote by ?/„ = x„(0) the starting position of particle n (non-random). 

We always start with a finite number of particles, but since the interac- 
tions are due only to the blocking from the right, it is effectively equivalent 
to have the index n varying over N. Limiting cases when n varies over Z can 
also be treated as appropriate limits of finite systems. For particle-dependent 
jump rates, we derived in [4] the general formula for the joint distribution of 
any subsets of particles at time t. This result is restated as Proposition HI 

In order to apply this result, we need to set the initial positions y^s and 
the jump rates v^s. In this paper we consider the first M particles to have 
jump rate a and the rest having unit jump rate: 

» = 2(M-A (2.1) 

The choice of setting the last a-particle at the origin is due to a simplification 
in the specific situation where we will take the M — > cxd limit. 

Macroscopic description for < a < 1 

The fluctuation results will depend on the macroscopic behavior, so let us 
first describe it. By macroscopic scale we mean when spatial directions are 
linearly scaled with time t. On that scale, for a G (0, 1), the effect of finitely 
many slow particle (s) is equivalent to having a starting density of rate 1 
particles on N equal to 1 — a. A particularly important case is M = 1, 
for which by Burke's Theorem [11], the initial condition is exactly equal to 
alternating deterministic on Z_ and Bernoulli- (1 — a) on N. 
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Let g{^, r) be the macroscopic density of particles, 

g{C,, t) = lim P(there is a particle at [C,t] at time rt) (2.2) 



t—tca 



Then, the average current from TASEP dynamics through position [^t] at 
time rt is given by ^(1 — g), from which it follows that g satisfies the Burgers' 
equation [29] 

drg + d^{g{l-g)) = 0. (2.3) 

To get the large-time macroscopic density one has to solve (I2.3p with initial 
condition 

The solution at r = 1 is as follows. For a G [0, 1/2), it has a discontinuity 
at ,^ = a — 1/2 (in this case, one needs to use a conservation law to obtain 
this solution, see e.g. [39]), 

.«.l) = |f' "MS- (2.5) 

^ I 1-a, if ^> (a -1/2), ^ ^ 

while for a e [1/2, 1] 

1/2, ife<0, 
g{^,l)={ (1-0/2, if^e [0,2a-l], (2.6) 
1-a, if^>2a-l. 

So, for large t, the density of particles in the lattice-scale, 

p{x,t) = P(there is a particle at x at time t) = g{x/t, 1), (2.7) 

see Figure [T] for an illustration. 

A consequence is that the number of particles moving with average speed 
a is around (1 — a)t/2 for a < 1/2 and (1 — a^t for 1/2 < a < 1. Also, the 
macroscopic position at time t of particle n = [ut] is 

a - z//(l - a), if z/ e (0, mm{^, (1 - a)^}), 
x,(z/) := lim t-^E(x[,t](t)) = { 1-2^, if z/ G ((1 - «)^ \), 



t— >oo 



1/2 - 2iy, if z/ > max{if^, i}. 



(2i 

where the second case occurs only for a G (1/2, 1]. 
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t = 



1 — a 



t > 1 1 



p{x,t) 



1/2 



t = 



X 



1 — a -6- 



t> 1 1 



(a - l/2)t at 



p{x,t) 
1/2 

p(x,t) 



H — ^ 

(2a-l)t at X 



Figure 1: Particles' density for large t: (a) for < a < 1/2, (6) for 1/2 < 
a < 1. The big dot is the position of the right-most slow particle. 



Process diagram 

As we have seen, there are different types of macroscopic behavior for the 
density. For example, when a G (1/2,1), there are two plateaux in the 
density joined by a linearly decreasing part. The plateaux are of different 
nature, since only the right one is influenced by the a-particles. So, the limit 
process of particles' position varies depending on which part of the process 
diagram the parameter are in, see Figure [21 

For keeping the presentation of the results as simple as possible, the 
results stated in the remainder of this section are the particularization to 
fixed time. However, the results hold in greater generality and span from the 
fixed time to tagged particle problem. The general statements are contained 
in the following sections. 

(1) Dyson's Brownian Motion region. For fixed time t let us consider particles 
with number n = [ut] with u G (0,min{i^, (1 — a)^}), i.e., we are in the 
right plateau with density 1 — a. In the diffusion scaling limit, the Mth a- 
particle has GUE(M) distributed fluctuations. So, to get a non-trivial limit 
for the particles in the jammed region we have to look at fluctuations with 
respect to the macroscopic behavior on the t^^"^ scale. Therefore, we set the 
rescaled process as 



X[M+ut] (t) 



[v)t 



-a{u)ty^ 



a\u) = a{l-v/{l-af), 



(2.9) 
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Figure 2: Process diagram. The black thick hne is the shock. On the dashed 
hne there is the Airy2--^i process. At the white dot there is the Airy2-^i,M,/t 
process. At the curved sohd hne the process is the AiryDBM^2- 

and Xq,(z/) = a + z//(1 — a), see (12.81) . Then 



DBM is the stationary process of eigenvalues of /3 = 2 Dyson's Brownian 
Motion on M X M Hermitian matrices, see Lemma [12] for a definition. The 
change in time — \n.a{v) is simply due to the non-stationarity of X((z/), while 
DBM is stationary. The complete statement is in Proposition [T3l 

(2) Shock region, M = 1. Consider now a G (0,1/2), where there is a 
macroscopic shock traveling to the left with speed (a — 1/2), and consider 
the important case of Bernoulli-(l — a) on Z^. as initial condition, that is 
M = 1. So, if a particle is already inside the jam, then it has t^/^ fluctuations 
with respect to the dashed line in Figure [31 On the other hand, particles can 
not move faster than they would in absence of the a-particles, in which case 
they fluctuate on a t^^^ scale around the dotted hne in Figure [31 So, on 
the t^/^ scale, the dotted line acts as a sharp cut-off and the result is the 



Proposition 1. Consider one slow particle with < a < 1/2 and the scaling 



limXt{u) = DBM(-lna(z/)). 



(2.10) 



following. 



n 



'-^t + r^t'l\ 
\t-2n- it 



1/2 



(2.11) 



For ^ > 




(2.12) 
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^1/2 




"''shock 



Figure 3: Illustration of the shock. The continuous line is the position of 
particles, the dashed line is the macroscopic position of particles inside the 
jam, while the dotted line would be the position without the a-particle. 
See [19] for an animation of the TASEP with and without a slow particle. 



where 



2 a(l — 2a) ^ 1 — 2a 



^ ^c=- V- (2.13) 

2(1 — a) 1 — a 

For^ < 0: 

lim F{xn{t) > a;(0) = 0. (2.14) 

Geometrically, ^ = are points on the dotted line, while C = are on 
the dashed line of Figure [3l This proposition is proved in Section 14. 2[ 

More precisely the picture is as follows: On the t^^^ scale, the random 
interface given by {(x„(t), n),n> 1} looks like a plot of a (non- homogeneous) 
Ornstein-Uhlenbeck process with average given by the dashed line. When the 
trajectory hits the dotted line, it sticks to it and does not fluctuate anymore. 
We identify this hitting position as the shock. This picture is consistent with 
Proposition [H A consequence of Proposition [1] is that the shock position 
has Gaussian fluctuations (for M = 1) with diffusion coefficient D = ■ 
This is the content of Proposition [161 The main novelty in this result is that 
we do not start with Bernoulli initial conditions on both Z_ and Z_|_ with 
two different densities and, more interestingly, we do not have to introduce 
second class particles to define the shock position, as it was the case for 
example in [12,16,17]. 

(3) Transitions and Airy processes. The last two new results related to the 
process diagram are the transition point (the white dot in Figure [2]) and the 
transition line between DBM and the Airy2 process (along the curved line in 
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Figure [2]) • The other remaining regions, are the ones where the influence of 
the a-particles is not present, so one clearly gets the same results obtained 
in [8]. 

The limit processes are different, but the scaling limit can be presented 
in the same way for all the cases. Indeed, consider n ^ ut with 



u > (1 -«)/2, for a G (0,1/2), 
z/>(l-a)2, for a G [1/2, 1). 

Then, the rescaled process for fixed time t is given by 



(2.15) 



XlM+ut-2rtV^]it) - x«(z/ - 2rt ^/^)t 
Xt{T) := . (2.16) 

The new transition process is at the macroscopic point given by z/ = 1/4 
and with a = |(1 + Kt~^^^). In Theorem [19] of Section [5] we prove that 

lim Xt{r) = S,A2-^i,M,.sAr/Sh) (2.17) 

i— »oo 

with = 2"^/^, Sh = 2"^/^, and r i-^ A2^i,m,k{t) is given in Definition \TE[ 
The second transition process is at the line v = {1 — a)^, for a G (1/2, 1). 
In Theorem [22] of Section [5] we prove that 

lim Xt{r) = S^Aj^BM^rlSh) (2.18) 

1 /3 

{l-a)a \ ' C _ (1-")^ c2 



with = (^ (i_a)(2-a) ) ^ = ^-^Sf„ and r t-^ ^dbm^2(t) is given in 
Definition [21] The process ^dbm^2(t) has appeared before, see [1,2,10,23] 
(with sometimes the time direction inverted). 

Finally, to complete the process diagram we state what the limiting pro- 
cesses are in the a-independent cases in Section [6] where one has either the 
Airyi, the Airy2 or the Airy2-^i process. The fixed time t results were already 
contained in [5]. 



TASEP with a reflecting wall 

The last result of this paper is of a different nature, since we do not have 
slow particles. Instead, consider the case a = 2 but with M = oo. Then 
the particle that started at the origin moves with average speed 1 and has 
fluctuations of order t^/'^. The other "normal" particles, have jump rate 1 and 
are blocked by the last a-particle, but in contrast to the jam situation they 
do not have the tendency of filling up the gap rapidly, since the "wall" moves 



9 



with their natural speed. Let us call particle n the one starting at —2n. 
Particle 1 fluctuates on a t^/^ scale, that means that from its perspective 
the last a-particle is like a moving blocking wall. Viewed from the "wall", 
particle 1 does essentially a reflected random walk in continuous time, and 
particle 2 does a random walk reflected on particle 1 and so on. 

Consider a sequence of particle numbers (not rescaled with time) and 
times ti = Tit. Since particles at time tj are approximately at position ti 
(speed one), we define the rescaled random variables 

t^Mt) = ^^^^^^. (2.19) 



We can compute the correlation functions of Xt{i)^s only if they are space- 
like, property denoted by ~ and defined by 

(ni,ti) ~ (712,^2) (ni,ti) -< (^2,^2) or (^2,^2) -< (ni,ti) (2.20) 

with 

{ni,ti) -< (^2,^2) ni < n2,ti > ^2, and are not identical. (2.21) 

Then, our result proven in Section [7] is the following. 

Theorem 2. For any given m = 1,2, . . ., let us choose m space-like couples 
{rii, Ti), 1 < i < m. Let pf^\^i, ■ ■ ■ , ^m) be the m-point correlation functions 
o/Xt(l),...,X,(m). Then 

^l_imp^)(ei,-.-,U =det [K-''''''{{n„ei),i,;{n„e,),Q]^^^^^^^ (2.22) 
where 6j = In(rj). 

The kernel i^^^iUE extension of the antisymmetric GUE minor ker- 
nel [22] defined as follows (see Lemma [2^ for an integral representation). 

Definition 3. The extended kernel K^'^^^ is defined by 

i^^«u^((ni,ei),6;(r^2,^^2),6) 

_ 2 . ^ sign(£)e-(^^-^^)^ 

where sign(£) := I if i > I, sign(£) := —1 if I < 0, and the interval of 
summation I is 

/ = {l,2...,L(n2 + l)/2j}, zfin,,9^)y^in2,e2), 

/ = {-oo,...,-l,0}, ^fin,,e,)^ {712,62). 

The functions Hk{x) are the standard Hermite polynomials with normaliza- 
tion J^dxHk{x)Hj{x)e~^ = 5k,jk\2'^ ^/t: . 
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In the fixed-n specialization of the i^^^UE j^g^^^el, i.e., for = n2 = n, 
this kernel becomes the one of a system of N non-intersecting Brownian 
motions (rescaled to become stationary) as follows: 

(a) if n = 2N — 1, the Brownian motions have the reflecting wall at the origin, 

(b) if n = 2A'", the Brownian motions have the absorbing wall at the origin. 

These processes were introduced and studied in [25,26] but the kernels were 
not exphcitly provided. Also, the kernel for n — 2N was obtained in the 
study of N Brownian excursions [35]. 

3 Determinantal structure and kernels 

We start by stating the general formula and then particularize to our choice of 
jump rates. Consider particles numbered by 1,2,.. ., with particle j starting 
from site yj and jumping to the right with hopping rate vj. The joint distribu- 
tions of particle positions are obtained as a specialization a{t) — t, b{t) — 
of Proposition 3.1 in [4]. To state the result, consider the set of numbers 
{vi, . . . , Vn} and let {ui < U2 < . . . < Uy} be their different values, with ak 
being the multiplicity of Uk- Then we define a space of functions in a;, 

K = span{x'<, 1 < A; < i/, < Z < ctfe - 1}. (3.1) 

The next statement holds for finite sequences of (distinct) events in the (n, t) 
variables which are space-like. 

Proposition 4. Let us consider particles starting from yi > y2 > ... 
and denote Xj{t) the position of jth particle at time t. Take a sequence 
of particles and times which are space-like, i.e., a sequence of m couples 
S = {(^T-fc, tfc), A; = 1, . . . , m I (nfc, tk) -< (rife+i, ife+i)}. The joint distribution 
of their positions {tk) is given by 

m 

P( n {^nuih) > Ofc}) = - XaKXa)e^{{{ni,h),...,{nm,tm)}xZ) (3-2) 

k=l 

where Xa{{'>T'k,'tk),x) — ^{^ < ^fe)- Here K is the kernel with entries 
K{{ni,t^),x^;{n2,t2),X2) = -<p^^''''''^'^''''''^\x^,X2)+K 

(3.3) 
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where 

k=l 

^ ' ^ 27ri Jr^ W ^/;^i+ni-=^2-n2 _ . . . _ 

(3.4) 

With V we mean {vm+i, • • • , ^712}- The contour Fq,^ is any anticlockwise 
oriented loop that includes and the elements of v. The functions '^^-j' 
j > 1 are given by 

The functions {^^-j}i<j<n o-f^ characterized by the two conditions: 

(K'-,, K'-k) ■= E Kt,{x)KU^) = kk, l<J,k<n, (3.6) 

and span{<l>^;^^(a;), 1 < j < n} — Vn- 

The notation 2^ dzf{z) here and below means that the integration 
path Fk goes around the poles of f{z) which are in the set K. 

In our situation, the orthogonalization gives the following result. 

Lemma 5. For n < M, the n orthogonal functions are 

1 / dw {w{w - a))"-^e*"' 

^"•^ (3.7) 



= 2^ ^%K.+i)((. + l)(. + !-«))-.>! + 2 - a), 

where j = 1, . . . ,n. For n > M + 1, we have two cases: 
(a)forj = M + l,...,n, 

f3 8) 
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(b)forj = l,...,M, 



X = i dv i dz—, 3.9 

(1 + z)^'+2»-2M ^ 

^ e*(-+i)(z(l + ((t; + + 1 - 

Our original derivation of the orthogonal functions was based on the 
known orthogonal functions for the case where all the particles have the 
same jump rate (see [6,30]), and then by employing Gram-Schmidt orthogo- 
nalization procedure. A similar procedure could be used also for more than 
one jump rate different from one, but we did not do it. However, once the 
orthogonal functions are determined, it is easier to verify the orthogonality 
by direct computation, and this is what we do below. 

Proof of Lemma\^ The formulas for \E'^^^(a;) is just a simple substitution of 
into dSISD. Notice that for x < 2M - 2n, = 0. Therefore, 

xeZ x>2M-2n 

in which x-dependent terms under the integrals are given by 

/ ^ _^ ^\x+2n-2M ^ 

y ( ) = provided \w\ > + (3.11) 

^-^ ^ w I w — [1 + v) 



x>2M-2n 



With these preparations we can prove the orthogonal relation needed by the 
theorem. 

Case n < M: We have 

If,! , iwiw - a)Y-^e^'"{2v + 2- a) 1 
dv (h dw- 



(27ri)2 /r_, /ro,+. {{v + l){v + I - a)Y~i+^e^^-+^) w - {v + I) 
--— i dv{2v + 2- a){{v + l){v + 1 - a)y-^-^ 

— / dzz^-''-^ = 5k,j, 
2vri iro 

(3.12) 
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where we used the fact that after summing over x the pole at w = disap- 
peared (since n—j > 0) and then the change of variable z — {v + l){v + l — a) 
with dz — {2v + 2 — a)dv. 

Case n > M+1: We have to do four computations, depending on whether 
j, k are larger or smaller than M. 

Case j,k > M + 1: We get {^n-j^ ^n-fe) ~ ^k,j by the same computation 
as in the n < M case but with a replaced by 1. 

Case j, k < M: Also in this case, after summing over x the pole Bit w — 
disappears but instead there is a simple pole at w — z+1. This can be easily 
integrated out and we get 

. ^"'*,\ 

(2vri)2/r._. /r„,„ {{v + l){v + I - 0))^'^+^ {z-v){z + v + l) 
= / d^;((^; + 1)(^; + 1 - a)y-''-\2v + 2 - «) = (5,-^ 

(3.13) 

where we used that after integrating out the w — z + 1 pole, the variable z 
does not have a pole at 2; = anymore. 

Case j > M + 1, k < M: In this case, the sum over x and then the residue 
at w = f + 1 leads to 

(3.14) 

because there is no pole at v — anymore. 

Case j < M, k > M + 1: It is slightly more tricky to check the orthogonal- 
ization in this case. After the sum over x and the residue at w = 2; + 1, we 
get 

\^n-j' ^ n-kl 

1 f 2v + 2-a r {z{z + l))^-''{l + 2z) 

^{2niyjr^_, \iv + l)iv + l-a)r-j+^fr,^ ' iz-v){z + v + l) ' 

(3.15) 

This time, both poles at ^ = and z = v contribute. One notices that the 
integrand in z has four poles in the whole complex plane: z — —1, 0, —1—v, v. 
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By the change of variable z = —1 — w, we have the identity 

, (z(z + l))^^-Hl + 2z) r , (w(w + l)y^-''(l + 2w) 

dz—, -r: ^ -— — = (h dw- 



ro,„ iz-v){z + v + l) /r_i,_i_, {w - v){w + v + 1) 

1 / ^ {z{z + l))^^-''{l + 2z 

dz 



2 /r_,,o._i_„.„ {z-v){z + v + l) 

(3.16) 

Since M — k < —1, the integrand is 0{l/z^) at 2; ^ 00, thus the integral 
fl3.16p is zero, which implies then {^n-j^ ^n-k) ~ ^- '-' 

For our analysis we will not focus on the first M — 1 particles, since it 
corresponds (up to a time-change) to the cases analyzed in previous works. 
Here we'll focus only on particles' positions of the ones with jump rate 1. This 
is the reason why in what follows we write the kernel only for ni,n2 > M. 
With the expressions of Lemma [5] we can rewrite the kernel (13.31) in the 
following way. 

Proposition 6. For rii, n2 > M +1, the kernel has the following expression 

K{{n,,t^),x,;{n2,t2),X2) = -0«"^'*^)'("^'*^»(a;i, X2) 

+ ir«((ni, ti), xi; {712, t2),X2) + K^^\{n^,t,),x,; (^2, ^2), ^2) (3.17) 

where 

(P ^ ^' ^'~2'Ki%^W y;^l+2ni-X2-2n2 H{n,M)<n2M)h 

(3.18) 

and 

7>(\\n N , ^ N ^ I . I dwe*i"'(w(w- 1))"!-^ 
H(ni,ti),xi;(n2,t2),X2) = j^dv — ^,^^2n,-2M 

eHv+^){v{v + l)Y^'^^ [w + v){w -V -I) 



and 



k^'^\{ni,ti),Xi]{ri2,t2),X2) = / dv l dz I — 

e*l"(w(w - 1)Y^-^^{W{W - a)f^ (1 + ^)^2+2n2-2M 

1 {I + 2z){2v + 2 - a) 

{{y + l)(t; + 1 - a))*^ (^ - i;)(2 + + l)(w - 1 - f )(w + 1 - a + t;) ' 

(3.20) 
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This proposition will be used in Section 

Proof of Proposition\^ For all > M, we have w^i+i = 1, thus fl3.4p implies 
that 0(("i'*i)'("2,t2))(^^^a,2) in ([33D equals Xa) plus the pole 

at f = 1 (we return to it shortly). For the rest, we divide the sum over k 
in (13.41) into the sum over [1, . . . , M] and the sum over [M + 1, . . . , na]. We 
define 



n2 



K^^\{n,Mlxr;{n,MU,)= ^ K%{^i)Kt\^^2) (3.21) 



k=M+l 

and 

M 



K(2)((ni,tO,Xi;(n2,t2),X2) = $^vI;-:^^^(:,,)$-:^^^(:,,). (3.22) 

A;=l 

Remark that ^^^%{x) = for k > n2 + 1. Therefore we can extend 
the sum in (I3.2ip to infinity. Then, if we take v small enough and w large 
enough, satisfying |f (w + 1)| < \w{w — we can take the sum inside the 
integrals. Explicitly, we get 

rWiw. ^ . N N ^ I . I dw; (w(?i; - l))"i-"2e*i'" 
H(ni,ti),xi; (n2,t2),X2) = f^^ y,^.+2n,-2M 

k=M+l ^ ^ " 

where the integration contours have to satisfy |f (f + 1)| < \w{w— 1)|. Then, 
using 

- (^(^ - 1))"2-^ _ (^(^ _ l)Y2-M 1 

f^i^^^{v{v + l)Y^-^+^ {v{v + l)Y^-^^ {w-{l + v)){w + v) ^' ' 

we obtain (13.191) plus the pole coming from w = v + 1. This contribution 
cancels exactly with the contribution of the pole at = 1 of 0. 

We now show that i^*^^^ = K'^'^\ For the computation of K^'^\ remark that 
the formula for ^^^-fcl-^) used for k < M gives exactly zero for k > M + 1. 
The reason is that the pole at f = a — 1 disappears. Therefore we can use 
the integral representations for k < M and extend the sum to = cxd. Then, 
provided that \{v + l){v + 1 — a)\ < \w{w — a)\, we can exchange the sum 
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and the integral, which gives 

rW2W. N . ^ ^ ^ I dw r r ^ («;(«;- l))"i-^e*i"' 

;i + 2)^2+2n2-2J\/ [i + 2z){2v + 2-a) 



X 



e*2(^+i)(z(;z + 1))«2-JV^ (2 _ + ^ + 1) 



00 



Then we substitute 

{w{w - a))^'-'' {w{w-a))^ 



((f + l)(t; + 1 - {{v + l){v + l- a))^^ {w - 1 - v){w + v + 1 - a) 

(3.26) 

into (13:251) to get ([320]). The condition \{v + l){v + 1 - a)| < \w{w - a)| 
is satisfied for any w G Tq^^ if we choose v sufficiently close to a — 1. Both 
poles w = a — 1 — V and w = 1 + t> lie inside Fo^q. Thus, K^"^^ is given by 
K^'^^ but with the poles for w = 0,a — l— v,l + v. Consider the contribution 
coming from the pole at w = f + 1, which is a simple residue. Computing 
this residue one immediately sees that the pole at w = a — 1 is not present 
anymore, thus the integral is zero. □ 

In the applications we'll use some special cases of the kernel too. In 
particular for M = 1 and M = 00 (with rii — M finite). Let us write the 
kernel explicitly in these cases. 

3.1 Special case: M = 1 

Corollary 7. For M = 1, the kernel has the following expression. For any 
ni,n2 > 1, 

K{{nuti),x,;{n2,t2),X2) = -0«"^'*^)'("^'*^«(xi, X2) 
+ ^ i dv ^ dwe*i'"(u;-l)"i-^(l + w)^2+n2-i {l + 2v) 



(27ri)2 Jr^ Jp^^ _^ w ^xi+m-i (,t2{v+i)yn2-i + _ ^ _ 1) 

1 / , e*i"'(w - 1)"!-^ 1 / , (l + f)^2+n2-i l + 2v 

H : (h dw ^ ■ — : (h dv 



2711 27ri7ro,<._i e*2('^+i)i;'^2-i (f + 1 - a)(f + a) ' 

(3.27) 

with (f) as in Proposition 0. 
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Proof of Corollary One simply substitutes for M = 1 in the expression of 
Proposition [HI Then, the integral over v around a — 1 is computed easily 
since it is a simple pole. □ 

The product structure of K^'^^ (the last term in fl3.27p ) is straightforward 
if one looks back at its definition fl3.22p . So, for M = 1 we have a rank-one 
perturbation. 



3.2 Special case: M = oo 

We want to get the M oo limit but with rii — M finite. Therefore, consider 
the kernel K{{M + ni, ti), a;i; (M + r;,2, ^2), 2^2) and take the M 00 limit. 
We define 

Koo{{ni,ti),Xi]{n2,t2),X2) ■■= lim K((M + ni, ti), Xi; (M + n2, ^2), 2:2), 

(3.28) 

and the limit kernel K^o is given as follows. 
Corollary 8. Forni,n2 > 1, we have 

K^{{n,,h),x,-{n2,t2),X2) = -0«"^'*^)'("^'*^»(xi, X2) 

Iff, e*i'"(u; - (1 + t;)^2+n2 i ^ 2v 

dv (b avj- 



(27ri)2 _^ yjx,+m+i (,t2{v+i)yn2 + _ ^ _ 1) 

-1 / , / , e*i"'(w; - (1 + '?;)^2+n2 i ^ 2v 

dvu (h dv- 



(27ri)2 Jp^ ^xi+ni+i ei2(i'+i)yn2 (^y I ^ (y){w - V - a) 

(3.29) 

wzi/i as in Proposition{^ 

Proof of CorollarylBi The only not straightforward term is K^'^\ For M ^ 00 
the pole at w = disappears and one just integrates out the simple pole at 
vu = a — 1 — V. The result does not depend on M anymore, namely 

-1 f , f , (l + ^)^2+n2gti(a-l-^;)(;^_2_^;)"l 1 + 2z 

dt> d) dz- 



(27ri)2 Jp^ Jp^ ^ e*2(^+i)z"2 (a - 1 - t;)^i+"i+i (^ _ ^;)(^ + ^ + i) ' 

(3.30) 

Changing the variable w = a — 1 — v and then renaming z with f we get the 
result of the statement. □ 

Notice that for a = 1 the combination of the two integrals in (I3.29p is just 
the residue at w = —v, which is the kernel for alternating initial conditions 
already obtained in [4,6]. Moreover, for ni = n2 = n, the kernel can also 
be seen as rank-ra perturbation of the kernel without the K^^^ contribution. 
This kernel will be used explicitly in Section [71 
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3.3 Modified kernel useful for the shock region 

For the asymptotic analysis in the case of finite M and a < 1/2, the shock 
situation, there is an interval around the shock for which the kernel has a 
diverging part in the t — >■ oo limit. This, however, does not mean that the 
system is ill-defined, because the distribution of the particles' positions is 
given by the Fredholm determinant of the kernel, not by the kernel itself. In- 
deed, we obtained a new kernel -ft'shock such that the Fredholm determinants 
agree. More importantly, in the new kernel (which is not a trivial conjuga- 
tion of K) the divergence disappears. To define A'shock let us introduce the 
following function. Set 



^ 1 {l + 2z)i2v + 2-a) 

{{v + l){v + 1 - a))^'-^+^ {z~v){z + v + l) ' ^ • ^ 

where w is either 0, —1, v. Then, for j = 1, . . . , M, we have 
Define the new kernel -ft'shock for ni,n2 > M as follows: 

+ K'^'^\{ni,ti),xi; {n2,t2),X2) + K^y^^^^{{ni,ti),Xi; (^2,^2), 3:2) (3.33) 

C2) 

where -K'shock defined as follows: 

^shLklK.i^O.a;!; (n2,t2),a;2) 

M M 
fc=l k=l 

(3.34) 

This means that for (ni, ti) -< {n2, t2) instead of the poles at 2; = 0, f in (13. 9p 
we have the poles at z = —1, f , and otherwise only the pole at 2; = f . The 
same changes in the poles will then occur in the triple integral representa- 
tion fl3.20p . This kernel will be useful for a < 1/2 because of the following 
Proposition. 

Proposition 9. For given Oi, . . . , G Z, we have 

(3.35) 

where Xa{{nk,tk),x) = l{x < ak). 
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Thus, by ( I3.2p the joint distribution of particles' positions can be com- 
puted with the kernel i^'shock instead of K. 

To prove Proposition [9] we'll use the following relations. 

Lemma 10. Forni,n2 > M. Then 

QZ~),-M)K^'\i^i^t^), xi; (^2, h),X2) = QZ%i^2), (3.36) 

Ql^'^\^,{x,)K^^\{n,,t,),x,- (^2, h),X2) = 0, (3.37) 

XI 

E^n::^^,,-i(^i)0«"^'*^^'^"^'*^^n:^i,^2) = QZ-U(^^)- (3-38) 

XI 

Proof of LemmaHE Recall the definition of fsT^^) and K'-^\ see fETD - flX^ : 
with /i = [M + 1, . . . , 77-2] and I2 = [I, . . . , M]. Then, we have to compute 

{QZ'-U^Kl'-\) = YQZ-Ui^)K['-U^), (3.40) 

so some of the computations are very close to the ones we made for the 
orthogonalization in Lemma 

Let us start with (13.371) . i.e., 1 < k < M. Then, the expression of Q is 
like $ in (13.91) but with Fq,^ replaced by r_i. Thus, compare with (I3.13p . we 

get 

Wni-j,-l' ^jii-fc/ 

1 / , / ((z + 1)(^ + 1 - a))''-" (1 + 2^)(2t; + 2 - a) ^ 
(27ri)2 7r^_^ Vr_, ' {{v + l)iv + 1 - a))'^ + ^ + 

(3.41) 

because the pole at 2; = — 1 vanishes for k < M. This implies (I3.37p . 

Next we prove (13.360 . i.e., k > M + 1. It is similar as before but with the 
^ taken from (13.80 instead of (13. 9p . Then, compare with (I3.15p . we have 

Wni-j,-l' ^ni~k/ 

(2vri)2 7r_, /r_, (^ - t;)(z + + 1) ((t; + 1)(^; + 1 - a))AW 

(3.42) 
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Thus, LHS of fl3.36p is given by 

"2 

E ^Qnt\-i^K%)Kt\{^2). (3.43) 

k=M+l 

We use the fact that ^^t.k{^) = if /c > n to extend the sum to infinity. 
Then, provided \w{w + 1)| < \z{z + 1)| we can take the sum inside the 
integrals; exphcitly we get 

„ , , 1 / ,1,1. (l + 2w;)(l + w;)^2+n2-Af 

(Km = 7 4> dv i dz (p dw- -f^-r '— 

^' (27ri)3/r_^ /p_^ ^ e*^(-+i)«;«2"A^ 

(l + 2z) 2i; + 2-a {w{w + 



X 



_ ^)(^ + ^ + 1) ((^ + l)(^ + 1 - a))*W ^Z^^^ + l))k-M 

, ^ , :i+w)^'2+"2-^^ 1 

df (b dw- 



(27ri)2/p^_^ Jp^^ e*2('"+i)'u;"2-^^ {{v + l){v + 1 - a)Y^-^+^ 

{l + 2w){2v + 2-a) 



{w — v){w + V + 1) 



(3.44) 



where we integrated the pole at 2; = w arising from the sum over k. This is 
however nothing else than Qnl'-jfli^'^)- 

Finally, we need to verify (13.381) . One divides the sum over xi in [0, 1, . . .) 
and (. . . , —2, —1]. Then use 

1 + z\ w 



xi>0 



= if l + ^<w 

\ w / w — 1 — z 



xi<0 



(3.45) 



= if\l + z\>\w 

\ w J w — I — z 



The two sums can be taken inside the integrals provided the contours satisfy 
once \w\ > |1 + -2| and the other time < |1 + zj. The integrands are the 
same up to a sign, which means that the net result of the sum is just the 
residue at w = 1 + z. Then (I3.38P easily follows. □ 

With this result we can now proceed to the proof of Proposition [HI 

Proof of Propositionl^ In this proof we let n stand for a pair {n,t) for short 
and write like i^nj^^jl^i? ^2) to represent ((ni,ti), (^2,^2)) block of K. Let 
us define 

M 
k=l 
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Since ^^tki^) = for x < 2M — n — k — 1 and Q^'lk ^liv) = for y > M — n, 

it follows Ejll ^n-k{^l)Qn\~l{^2) = Q li Xi < M - U - I OI X2 > M - U. 

Hence 5'„,„(xi,X2) = for X2 > xi, i.e., 5* is lower triangular with diagonal 
being zero. Hence to prove the proposition, it is enough to show 

[1- S){1- K) = 1- k. (3.47) 

Since 5" is block-diagonal, the ni,n2 block of fl3.47p writes 

(1 — >S')„^,„^(1 — K)ni,n2 = (1 — -f^shock)ni,n2, (3.48) 

which is proven using the result of Lemma [TU] as follows. We use the notation 
K = K^^^ + below. 

Case rii = Then LHS of (13.481) is 

(1 *S'jji,ni)(l -^ni,ni) 1 Sn-^^^m -^7ii,ni ~l~ 'S'^j^ ,ni -^ni ,ni • (3.49) 

We have Kn^^m given by 

M 



fe=i 



and by Lemma [10 



M 



Sni,niKni,ni — ''^'^ ni~kQni-k,0- (3.51) 
k=l 

Putting together these relations we obtain exactly RHS of (13.481) . 

Case 111 7^ 77-2 and rii n2: LHS of (I3.48P is in this case given by 

(1 Sn2,n2) n2,ni -^712, ni ~l~ Sn2^n2-^ n2,n\- (3.52) 



Using 



and 



we get 



M 



Kn2,n, = + E Kt\(Qlt\o + Qlt\vl (3.53) 



fe=l 

M 



5'„2,n2-^n2,ni — ^n2-fcQn!-A;,0 (3.54) 



fc=l 



M 



[1 Sn2,n2)Kn2,nx — -^ia!™! "^^^ n2-kQ ni-k,v (3.55) 



k=l 
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which is the claimed result. 

Case rii ^ n2 and rii -< in this case, LHS of f l3.48p is given by 

(-W — Q V— — _|_^("l."2)_|_ C '17 _C i(jll,n2) 

y-"- >-'ni,niJ\ ''■^ni,n2~'r J -^^ni,n2~'r ~>-''ni,ni-^^ ni,n2 >-'ni,ni't' 

(3.56) 

This time Lemma [TO] tell us that 



M 



^^ni,n2 ^^ni,n2 ^ / ^ ^ ni-k\^ n2-k,0 ^ ^n2-k,vJi 
k=l 

M 

Sni,niKni,n2 = ^ ^ni-fe^na-fc.O' (3.57) 
k=l 
M 

fc=l 

These relations imply the claimed result. □ 
3.4 Special case: M = 1 

For later use we explicitly state a corollary of Proposition [9l For M = 1, 
the extended kernel K is given in Corollary [71 Proposition [9] tell us that 
the Fredholm determinant can be also computed using the modified kernel 
-f^shock, which has the same expression as fl3.27p but with the last term the 
integration for v is around the poles —1, a — 1 instead of 0, a — 1. In particular, 
for (ni,ti) = (^2,^2) = {n,t), the kernel is given as follows. 

Corollary 11. For M = 1, the one-point modified kernel is given by 

^shock((n, t), x; (n, t),y) = i^„,t(x, y) + f{x)g{y), (3.58) 



where 



Kn^t{x,y) = ——r (h dv 



d«;e*"'(w; - 1)"-^ (1 + f)2'+"-^ {l + 2v) 

(27ri)2 /p^ /p^ _^ W (,t{v+l)yn-l (^y + ^) (^y _ ^ _ 1 j 

1 / ^ e'^(w-ir-^ 
j[x) = -— (p aw- 



27ri /pjj ■njx+n 

If ^ (l + t;)2'+"-i l + 2t; 

9{y) = TT^ f dt; 



27ri7p_^^_^ gt(t«+i)^n-i + 1 _ ct)(t; + a)' 

(3.59) 

This result, together with Proposition [9] will be employed in proving 
Proposition [H 
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4 Jam regime 



Consider the semi-infinite system with 1 slow particle. By jam regime we 
mean the following two situations in which particles with jump rate 1 are 
slowed down by the slow particles: 

(1) for 1/2 < q; < 1: at large time t, the macroscopic density is continuous 
and has a plateau with density 1 — a. The plateau corresponds to the first 
{1 — a)H particles moving with speed a. 

(2) for < a < 1/2: in this case, the slow particle create a macroscopic 
shock and the density has a jump from 1/2 to 1 — a. Particles in the region 
of higher density move with speed a. The shock has a drift velocity equal to 
Vs = a — 1/2 (i.e., it moves to the left). 

With the results for case (2) we'll also be able to determine the law and the 
diffusion coefficient of the shock without introducing second-class particles. 



4.1 Fluctuations in the speed a region 

For large time t, particles with a particle number n < min{i^, (1 — Q;)^}t will 
move with the speed of the slow particles and will be very much correlated 
with the first M slow particles. What happens is that the Mth particle 
fluctuates according to the largest eigenvalue of DBM. Intuitively, then the 
other particles have jump rate 1, which is strictly larger than a, so that they 
fill the gaps more rapidly than if the jump rate would have been a. In doing 
so, their fiuctuation will be well correlated with the last slow particle, just 
shifted in time. Therefore one might expect to see DBM. 

Before stating the result, we define the limit object we'll get in the large 
time limit. The matrix- valued (stationary) Ornstein-Uhlenbeck process on 
M X M hermitian matrices also known as Dyson's Brownian Motion, DBM. 
It is the Markov process with transition density given by 

1 / Tr(M2-e-(^2-ri)j^^N^ 

Pin, M,; r„ M,) = ^.^^^ _ ^_.(.._.,)))m ^ 2(1 - e-^(^^-n)) 

(4.1) 

for T2 > Ti, and the reference measure being flat over the independent 
entries of M, i.e., dM — J^^^lj^ dMj ^ (iReMjjdImMjj. Its finite- 

dimensional distributions are given by the following Predholm determinant. 

Lemma 12. For any given ri < . . . < G M, the joint distribution of the 
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largest eigenvalue of the stationary DBM process are given by 

P ||f|DBM(r,) < Sk^ = det(l - X.i^™X.)L^({n,...,r^}xM) (4.2) 
where the kernel is given by 

f (a;2 -a;ie-(^2-ri))2\ 
^^P 1 Km _2( _ -)\ I M-i 

(4.3) 

where pk{x) = Hk{x/ \/2)'K ^^'^2 '^/^(/c!) and Hk{x) is the standard Her- 
mite polynomial of degree k (see e.g. [27]). 

This result can be found in [24] with a shghtly different normalization, 
an extra \/2 in the space variable. 

Proposition 13. Let 7r{9) be a real-valued function on M with |7r'| < 1. 
Define the scaling 

t{e,T) = inid) + 9)T, 

nie,T) = M+[i7ri9)-9)T], ^ ' ^ 

where T is the large parameter. For < n < min{^^, (1 — i.e. for 

< 7i{9) < min{|^, ^a{2-af }^' '^^ ^'^^ inside the region of speed a. The 
rescaled process 



Xn,eTM9,T))-{at- ^^^'f^-^^ ) 

Xt{9) = ^ ' ^ " — ^ ^-^ — (4.5) 

—ayT 

where = a{TT{9) + 9) — ^^^z^P-- Then, in the large T limit Xt converges 
to the DBM process: 

lim Xt{9) = DBM(r(^)), with t{9) := -ln(a) (4.6) 

T—^00 

in the sense of finite dimensional distributions. 

Space-like paths include as particular cases: (a) fixed time with t = T is 
obtained setting 7r{9) — 1 — 9, and (b) fixed (tagged) particle with n — T by 
setting 7r{9) — 1 + 9. For more explanations about space-hke paths see [4]. 
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Path for Paths for K^^^ 



a ^+,2 a 




Paths for K^^\ part (a) Paths for K^^\ part (b) 

Figure 4: Steep descents for the different terms in the kernel. The path w 
passes either through a or at distance eT~^/^ from it. 



Remark 14. In the following proof, as well as in the others on asymptotic 
analysis, we present only the most important ingredients. First of all we 
state explicitly the steep descent path used for the analysis and the local 
series expansions around the critical points (from where the non-vanishing 
term arises). These two are the building blocks for the convergence of the 
kernel on bounded sets, for more details on the procedure see e.g. Lemma 6.1 
in [3]. We do not however prove convergence of the Fredholm determinants, 
for which bounds on moderate and large deviations are needed, for a simple 
example on how to proceed, see Lemma 6.2 in [3]. 

Proof of Proposition UR The result is obtained by analyzing the rescaled and 
conjugated kernel 

Kne,,^,;e2,^2) = ^P^a^T'/'Kiinuh),xv,in2,h),X2) (4.7) 
Conji 

with Hi := n{6i, T), ti := t{9i, T), and 

Xi := aU - 7^ - ii(TiT^'^. (4.8) 
1 — a 

The conjugation factor is given by Conj^ := e°'*'(a — 1)"'*/q;^'"'""'\ We use the 
kernel K in Proposition [6l 



26 



Let us start with (f). In this proof, define the notation = 7r(0j) — 9i and 
Ui = Ti{6i) + 6i. Then, with the above scahng, for -< (^2,^2), 

0«"i'*i)'("^'*^))(xi,X2) = — / Awe'3o{^^+'''''^.(^^+^2(^\ (4.9) 

27ri 



ro 



where 



/ a 

go(w) = (ui - U2){w - alniw)) + (ai - 02) ln(w - 1) + :; \n{w] 

\ 1 — a 

9i{w) = (^iCTi - 60-2) ln(w), 

^2(w) = -ln(w). 

(4.10) 

The condition {ni,ti) -< (^2,^2) means that 02 — Oi > and Ui — U2 > 
(at least one of the two inequahties being strict). The critical point for 
steep descent of go{w) is at w = a and as steep descent path we use 
Tq = {ae^y,y G (— 7r,7r]}. Indeed, 

^Re(^o(^ = ae'y)) = - (^{u, - u,) + «sin(i/) (4.11) 

which is negative for sin(?/) > and positive for sin(|/) < 0, so the maximum 
of Re((7o(u^)) is at w = a and the path Tq is steep descentl^l. By usual 
steep descent analysis, the relevant contribution of the integral comes from 
a ^-neighborhood oi w = a. There, we can apply Taylor series: 

g^iw) = g^ia) - i^{a\ - al) + ^(y^), 

g^{w) = gi{a) + ii/(6^i - 6^2) + 0(2/'), ^^'^^^ 
g2{w) = -ln(a) + 0(y). 

Notice that exp(T(yfo(a) + T^/^(7i(a)) = Conj]^/Conj2. Then, deleting the 
corrections 0{- ■ ■ ) accounts for an error of order (9(T~^/^) times the leading 
term. The change of variable y := zT'^/"^ implies that the leading term is 
given by 

— — / dzexp (-(0-^ - 0-^)2:72 + i2(6cri - 60-2)) 



Conji 11 f (6 - 6o"2/o"i)^ 

exp 



Conj2aiTV2^27r(l-a2Va?) "V 2(1 - aiK) 



(4.13) 



^For an integral I — dze^^''^\ we say that 7 is a steep descent path if (1) Re(/(z)) 
is maximum at some zq g 7: Re(/(z)) < Re(/(2;o)) for z e 7 \ {zo}, and (2) Re(/(z)) is 
monotone along 7 except at its maximum point zq and, if 7 is closed, at a point zi where 
the minimum of Re(/) is reached. 



27 



Now consider K^^\ With the above scahng, and after the change of 
variable v = z — 1, we get 

Iff (,T}o,i{^)+T^/^hAw) 2z-l (4.14) 

d^ (b dw- 



(27ri)2 7r^_^_^ e^/o,2W+Ti/2/i,2{.) + ^ _ i)(^ _ z)w 

with 

/ a 

fo,i{w) = Ui{w — a \n.{w)) + at I ln(w — 1) + ln(w^ 



I- a ' 7 ' (4.15) 



/i,j(w) = ^jdi ln(w;). 



Let us look for the critical points of /o,i and the steep descent paths. We 
have 

dw w(w — 1)(1 — a) 

that is, the critical points are 

uj^ = a and u;_|_j = l — -. (4-17) 

Ui[l-a) 

Using the relation 

< tti < Mimin{(l - a)/2, (1 - af} (4.18) 

one verifies the following relations: 

< = a < cj+i < 1, and > 1/2. (4.19) 

As steep descent path we choose: Vi = {z = 1 — re^'^,r = 1 — co'+_2} and 
To = {w = ae^'^}. Let us check the steep descent property. We have 

A Re(/„,(u>)) = -^^1^ (|„ - - (4.20) 

dip \w — 

which is decreasing while moving away from the critical point w = a, since 
by (14.181) the term in the parentheses is strictly positive. Also, 

— Re(-/o,2(^)) = rr^ [\z\ -a [I (4.21) 



dip ' \z\^ \ \ ^2(1 — a 
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which is decreasing while moving away from the critical point z = 
Indeed, using f l4.18p and fl4.19l) the term in the parentheses is strictly positive. 
Hence the 



Conj 



integral (14.141) is of order 
1 exp (TRe(/o,i(cj_) - /o,i("))) 



Conj 
Conj 



Conj 



12 exp (TRe(/o,2(^+,2) - /o,2(a))) 

!iexp(TRe(/o,2(a) - /o,2(^+,2))) = ^^exp(-5T), 6 > 0. 
l2 Conj 2 

(4.22) 

Indeed, /o,2(w+,2) > /o,2(a), because > a and -^fo,2iz) > on 

z e (a,a;+^2)- Therefore in the T — >■ oo limit, the contribution of K^^^ (con- 
jugated and rescaled) goes to zero. 

Finally, we need to consider K^"^^ of Proposition [6l We change the vari- 
ables V = V — 1 and z = z — I'm. (13.201) and get 

N / 1 f f f dw (w(w — a))^ 



(2vri)Vr. /r,, w {v{v - a)r 

{2z-l){2d-a) e^^°-^("')+^''"^M("') 

^ {z + d - l){w - v){z - d){w + V - a) e^/o.2(i)+ri/2/i,2(5) • ^ ■ ) 

We have the following two contributions: 

(a) : z around 1. From the above analysis on the steep descent paths, 
we choose Iw — a| = e/2 for e small enough, {w = (a — e)e^'^} and 
{z = 1 — re^^,r = 1 — uj^^2}- Since /o,2(i^+,2) > /o,2(«), we can choose e 
small enough such that the overall contribution is ^^^exp(— 5'T) for some 

6' > 0. Therefore in the T ^ oo limit, this contribution of K^'^^ (conjugated 
and rescaled) goes to zero. 

(b) : z around v. Integrating out the simple pole z = v we get a contribution 
equal to 



dv 



dw e^/o,iW+TV2/i,i(^) ^^^^ _ ^^^M 2v-a 



(27ri)2 /r^^_, w e^/o.2(s)+Ti/2/i.2(s) ({)({}- a))Af (^^ ^ y)(^^ + y ^ a)' 

(4.24) 

The integration paths are now chosen as Fq, = {|{; — q;| = i?T~^/^} and 
^o,a-v = {w = {a — LT~^/'^)e^y ,L > R}. With the above computations 
we have that To^a-a is a steep descent path, so that its leading contribution 
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comes from a T neighborhood of y = 0. There we can use Taylor series: 



2a2 

f f \ f ^ ^ , - a) .2^ (4.25) 

a 

iw{w - a))^'' = -a- 1)*'(1 + - «))). 

Denoting u = a(l + l^cxf^T-i/^) and w = a{l + Wa^^T-^/"^) we have that 
the leading term of (I4.24p for large T is given by 

^ dl^ / dVT- 



Conj2 fTiTV2 (27ri)2 /f^i^^ y^^ W -V e^'(-2/-i)V2+V6-2M 

(4.26) 

with L > R. By the change of variable W —W and —V we obtain 
(up to factors a/2 due to the different space-scaling) the extended Hermite 
kernel, see (2.13) of [24], which can then be rewritten in terms of Hermite 
polynomials. □ 



4.2 Fluctuations around the shock 

In this section we consider a < 1/2 and focus around the shock position. We 
consider the case of M = 1 slow particle, since it is also physically the more 
natural situation. Indeed, the system with one slow particle is equivalent to 
having stationary initial condition on Z+. The shock position, which at time 
t will be around position Vgt = (a — l/2)t, fluctuates on the t^^^-scale. How 
does it happen? The fluctuations of particles before entering the shock region 
live on a t^^^ scale, thus viewed from the t^^^-scale, these particles essentially 
do not fluctuate. So, one has some probability that the particle is not inside 
the shock region, which will then be very well localized, and if the particle 
is in the shock region, then it follows fluctuations of the slow particles. The 
occurrence of these two distinct intermediate scales allows us, in particular, 
to determine the diffusion coefficient of the shock. The result agrees with 
the argument in [33] modified appropriately for our situation, see below. 

The first result was stated in Proposition [1] and in order to prove it we 
recall the following result from [8] . 

Lemma 15. Consider the kernel without the slow particle, i.e., Kn,t defined 
in lis. 5(A) . and the rescaling 

n = z/t, Xi = \t-2n- C^^^/^ with v>l/A. (4.27) 
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Then, uniformly for Q in a bounded set, 

lim f/'Kr,,tixi,X2) = ir^,(Ci,C2), (4.28) 

where Kji,^ is the Airyi kernel [6, 18, 30]. 

With = we mean equivalent, since indeed to get a well-defined limit one 
has to do a conjugation of the kernel Kn^t- 
It is also known [21] that 

Fi(2s) = det(l - Kj^jLHis,oo),dx), (4.29) 
with Fi the GOE Tracy- Widom distribution function [34]. 

Proof of Proposition [11 Let us start with ^ > 0. From Proposition [9] and 
Corollary [TT] we have 

P{xn{t) > x) = det(l - XxKn,tXx - Xxf ■ gXx) 

= det(l - XxKn,tXx)il - (gXx, (1 - XxKn,tXx)~^Xxf)) 

(4.30) 

with Kn^t) f and g being defined in fl3.59l) . and Xx = l(-oo,a;). Compare the 
scaling flTTTj) and g2ZD: ^ ^ ^+r}t-^''^ and ^ ^t^'^- So, for any ^ > 
we will actually focus on the upper tail of the Airyi kernel and of the related 
distribution, i.e., Kn,t{x,y) ^ as t ^ oo (with conjugation), and 

lim det(l - XxKn,tXx) = 1. (4.31) 

t—^0O 

For the second term, we use 

(1 - XxKn,tXxy' = 1 + XxKrr,tXxil " X.^n.tXx)"'. (4.32) 

Thus 

{gXx, (1 - XxKn,tXxy^Xxf) = (gXx, f) + (gXx, (l - XxKn,tXxy^XxKn,tXxf) 

(4.33) 

and the fact that Kn,t implies then that the last term goes to zero. So 
lim P(a;„(t) > x) = lim 1 - {gXx, f) = lim {gXx, /), Xx = l[x,oo)- (4.34) 

t—^OO t— +00 t— >oo 

In the last step we used the orthogonality between g and /, namely {g, /) = 1. 
Under the scaling fl2.11l) . x(^) -|- n ~ ^at > 0, which means that the pole at 
V = —1 in the function g{y) defined in fl3.59p vanishes. Therefore, 

^x+n—l 
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and then 

rt-l 



dW(-)fi^y"V^)^^"^ (4-36) 



' |ui|>o 

2vri /|^|>„ w - a 



The steep descent path of F passes by the saddle point at w = a, but since 
it is a pole we have just to deform locally on a t"^/^ scale to pass on its right. 
The leading contribution is coming from a )f:~^/^-neighborhood of the w = a. 
Setting w = a + iyat'^^"^, we get 

F{w) = + iy{^ + a + 0{yh-''^). (4.37) 

The 0{y^t^^/'^) term is controlled by the quadratic term, and in the end we 
obtain 

\im{gxxj) = — dy , Ve > 

= -=L=/ ^dzexpHV(2a^))- 

Now consider ^ < 0. We first use a probabilistic argument. It is quite 
clear (by a simple coupling argument) that 

P{Xn{t) >X)< P{Xn{t) > X) (4.39) 

where P is the measure of the system without the slow particle. For the 
system without slow particle we have the result of Lemma [151 which tells us 
that 

P(x„(t) >x)= det(l - XxKn,tX.) ^0, t^oo. (4.40) 

The reason is that < corresponds to Q = ^t^^^ — oo as t ^ oo. Then 
f l4.40p follows from the non-degeneracy of the distribution Fi (no mass is lost 
at — oo). □ 

It is a bit more natural to look at the fluctuations with respect to the 
dashed line in Figure [31 Then, the result of Proposition [H rewrites as follows. 
Let 

F(0 := \im P(x„=[(i_,),/2+,ti/2] W >at-n- ^f/^). (4.41) 
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Then, F{^) has a jump at ,^ = namely 

1 1 r°° 

(4.42) 

This can be used to determine the diffusion coefficient of the shock without 
having to identify it with second class particles. When a < 1/2, the macro- 
scopic density has a jump from 1/2 to 1 — a. As we saw in Proposition [T3l 
before the shock the fluctuations becomes asymptotically Fi-distributed on a 
t^/^-scale, while inside the shock region are Gaussian on the t^^^-scale. Thus 
the position of the shock itself is localized on the t^^^-scale, see Figure [3] for 
an illustration. 

The question we want to address is how to determine its law and in 
particular its diffusion coefficient. Denote by Xshock(^) the position of the 
shock at time t. 

Proposition 16. In the large time limit, the shock is Gaussian distributed 
with diffusion coefficient D given by 

D = (4.43) 

1/2 — a 

In other words, 

1 

lim P(xshock(t) > (a-l/2)t-z/t^/2) ^ / da;exp(-a;V2/^). (4.44) 

Proof of Proposition {TR To prove the result we first have to understand what 
Proposition [1] says. Consider the particle with number n and look at position 
X rescaled as in 02.111) . The condition ^ > means that x is on the left of the 
dotted line of Figure [3] by — ^c)^^^^- Moreover, before reaching the shock, 
particles fluctuate only on a t^/^-scale away from the dotted line. Thus, for 
^ > ^c, Xnif) < X imphes that particle n already reached the shock. On the 
other hand, if particle n did not reach the shock region yet, then (on the 
t^/^ scale) it has to be on the dotted line (can not be farther to the right 
because of 02.141) ). Therefore, the probability that particle n has not yet 
reached the shock (i.e. a;shock(^) > Xn(t)) is equal to the mass at ,^ = ^c- 
Thus, from fH:i2l) it follows that 

lim P(xshock(t) > (a - l/2)t - z/t^/') (4.45) 



t—^OO 



2 poo 2 /"CXD 

= ^E=^ / dyexp(V/2^') = ^== / dxexp(-xV2D) 

after a change of variable. □ 
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Remark 17. The above argument is quite flexible and one could extend to 
the case of M slow particles instead of only one. We expect the following. 
Proposition [T] would be similar up to the distribution in (12.121) changed from 
Gaussian into the GUE(M) (the distribution of the largest eigenvalue of 
M X M CUE matrices) and the shock will have a GUE(M)-distribution with 
appropriate parameter, by the change of variable as in (I4.45p . 

This result can also be explained with an heuristic argument, following 
arguments in [33]. In the continuum hmit the particle density pt{x) is de- 
scribed by the viscous Burgers equation with noise (see (5.37) of [33]), 

Here e is the lattice constant, u is the diffusion constant and Jt{x) is the 
random current. The initial condition is divided into two parts, 

Poix) = psix) + Ve^ix). (4.47) 

Here ps is the deterministic part, 

/ \ i P- '■= 1/2, a; < 0, /, ,„^ 

Psix) = < ' n (4.48) 

^ \ p+ := 1 - a, a; > 0, ^ ^ 

and ^{x) takes into account the randomness in the initial conditions for x > 0, 

{^{x)^{x')) = I °' ... ^ °' (4.49) 

Note that in our present case there is no randomness for x < 0. The solution 
to (14.461) is of the form, 

pt{x) = ps{x - Vst - V7Db{t)) + O(v^), (4.50) 

with b{t) is the standard Brownian Motion. The shock front remains sharp 
but its center performs the Brownian Motion. The diffusion coefficient of the 
shock D is of our interest. 

Let us suppose that the initial density fluctuations move with constant 
velocity towards the shock and that this is the source of randomness of the 
shock location. At time t the particle density fluctuations starting from the 
region [vst,—Vst] have arrived at the shock so that ^'^ ^{x)dx represents 
the excess amount of particles comparing to the deterministic part. Since 
the difference of the density to the left and the right is p+ — p_, we would 
have 

^/Dh{t) = / ^{x)dx. (4.51) 

P+ " P- J Vst 
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Using fOUj) we get 



D = /!±il^ = (4.52) 
p+-p_ 1/2- a ^ ' 



This is the same as (14.431) . 

5 Transition processes 

In this section, we first focus around the critical parameter a = 1/2 and later 
on the Airy2 to DBM(M) transition. For a = 1/2, on a macroscopic scale 
the density is constant and equal to 1/2. However, the fluctuations to the 
left of the origin live on the t^^^ scale, while on the right they live on the 
scale. Here we consider a - 1/2 = 0{t~^/''^) and n - t/4 = 0(^2/3). We 
keep M fixed and finite. 

As before, we are not obliged to stay on a fixed time, but we can consider 
a space-like path described by a function 7r(^) with |7r'| < 1. Consider the 
space- like setting as in Proposition | 



t(r, T) = {n{e - tT-'/') + 6 - tT-'/') T, 
n(r, T) = M + [n{e - tT-^'^) -{6- rT~^/'^)] T, 

with > fixed ancj^ n{9) = 56/3. This ensures that macroscopically we 
focus at the transition region, which for a = 1/2 is around n = t/A. 

Here we consider a not necessarily exactly equal to 1/2. Instead, let us 
define 

1(1 + kT-'/'). (5.2) 



a 



Then, the rescaled process of particle position is given by 

' Xn{t) - {k - 2{n - M)) 



-Ti/3 



if n > t/4. 



_rpl/3 ' 11 ^ t/^. 



(5.3) 



In the large-T limit, Xt will converge to a well-defined limit process, 
•^2^1 M K, which we now define. 



^This does not mean that the function 9 i-^ Tr[9) is identically equal to 59/3, only that 
at 6* = 6* its value is equal to 59/3. 
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-72 7i 72 7i 72 

Figure 5: Illustration of the integration paths defining the kernel K^lf'^^. 



Definition 18. Let us set 



(5.4) 



The process A2^i,m,k is the process with m-point distributions at ti < T2 < 
. . . < Tm given by the Fredholm determinant, 



P n {^2^i,MA(^fc) < ^k} = det (1 - XsK^JXs) 



\k=l 



L2({ri,...,T,„}xII 



(5.5) 



where Xs{^k,x) = l{x > Sk)- The kernel is defined by 

K%Tir,, s,; r2, s^) = K'^^^<^\r,, s,; r^, ^2) + i^I^^^'Vi, ^1; ^2, ^2), (5.6) 
with 



ir*--«(n,Si;r2,S2) 
1 



^47r(r2 - ri) 



exp 



4(r2 - ri) 



l[r2>ri] 



+ 



(27ri)2 



□ Ulf /3+T2 - S2liJ2 



2W2 



dw2 / dwi — ^— 5-^ — - — 

^2 7^1 g»^f/3+Ti«;i-.i,«i (y;^ _ M;2)(wi + ^72) 



(5.7) 



and 



7-^trans(2) / \ 1 

Km,k {ruSi;T2,S2) 



(27ri)3 



gi02/3+T2)i;|— S2UI2 

dw / dw2 I dwi — ^— 5—: — 

' ' gUlJ/3+TlU)J-SlUll 



72 >^7i 
2W2 



M 



. (5.8) 



(W2 — U){W2 + u){Wi — u) \ U — K 

Here 72,72 : e'^'/^oo e~'^'/^oo, 71,71 : e^^'^^^'^oo and (7oes 

around only the pole at u = k anticlockwise. Moreover, —72 C 71, and 71,72 
passes on the left ofT^, (see FigurelEfor an illustration). 
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Theorem 19. The process Xj- defined in Ii5.3\) converges to the process 
A2->i,M,K, Tnore precisely 

lim Xt{t) = S,A2^^,M,nS. [r/Sh) (5.9) 

T — >oo 

in the sense of finite- dimensional distributions. The scaling coefficients Sh 
and Sv are given by 

49 \ 1/3 



In the fixed time case, t = T, we have 7^{0) = 1 — 6 and tt{6) = 56/3, 
from which 6 = 3/8, i.e., = 2'^/^ and St = 2'^/^ 

Remark 20. When M = we have K^^^^^'^^ = and the transition process 
is the Airy2^i, ^2^1, discovered in [8]: 

A^i,o,.(r) = ^2^i(r). (5.11) 

Proof of Theorem [Wi To prove the resuh, we have to analyze the large-T 
limit of the kernel in Proposition [6] under the following scaling: 

t. = fT-n{n'{6) + l)T'/\ 



n, = M+-T- nin'i6) - l)T'/^ (^-^2) 



2 

where 



Xi = \ti-2{n,- M)-%T^'^, 
Si, iiTi>0, 

s.-TfS,S-^\ ifr,<0. ^^-'^^ 

Higher order in the development of iiid — TiT"^^'^) are irrelevant since they 
corresponds to a T~^/^ perturbation of ^'{6). 

Then, we have to consider the rescaled and conjugated kernel 

i^r'(ri,si;r2,S2) := ^^^T^/''K{{n,,ti),x,- {112,12), x^) (5.14) 

with Conji := e*'/2(-l/2)"'-^^(l/2)"(^'+'^»-^). We need to show that 

hm Kr{r„sv,r2,S2) = S;;'K'^^:i{n/ S,, sj S,;t2/ S,,, S2/ S,). (5.15) 

The first two terms of the kernel f l3.17p are independent of a and their 
sum is the kernel without slow particles. This kernel was already analyzed 
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(a) (b) 

Figure 6: Steep descents for (a) K^^^ and (b) K^'^\ They satisfies n/A < q < 
q' < 7r/2, 7r/4 < g", q'" < tt/2 and the local modifications around the critical 
point at 1/2 are all only on the T^^/^ scale. 



in great detail in [8] with the slight difference that the space-like setting 
introduced in [4] was not known yet. However, at the level of asymptotic 
analysis there are no relevant changes. Thus here we just indicate the key 
steps. 

Let us first consider K'^^^ defined in (13.191) . After the change of variable 
V = V — 1 it writes 

1 f ^ r 2v-l e'^Mw)+T^/^fi,i{w)+T^/^f2,iM 

(5.16) 

with 

/!,,(«;) = -r,(7r'(^) + l)w - n{n'{e) - 1) \n{w - 1) + i(3 - n'{e)) \n{w), 
f2,M = Si\n{w). 

(5.17) 

From the analysis of Proposition 4 of [8] we have that the steep descent paths 
in (15.161) are chosen as illustrated in Figure El Next, the Taylor expansion 
around the double critical point of fo{w) , which is at w = 1/2, are given by 

foiw) = /o(l/2) - + o{{w - 1/2)% 

h,{w) = h,{l/2) - r,(5 - ?,iT'{d)){w - 1/2)2 ^ (^((^ _ 1/2)3), (5.18) 
f2Aw) = /2,.(l/2) + 2%{w - 1/2) + 0{{w - l/2f). 

The leading contribution to the kernel comes from the T~^/3-neighborhood of 
the critical point. The conjugation terms are just the value of the exponential 
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factor evaluated at the critical point. The 0{- ■ ■) term accounts into an error 
(9(T^^/'^) smaller than the leading one. Then, after change of variable 

w = i + ^WiT-'/yS,, v = l + \w2T-^'^/S,, (5.19) 

we get the final result 

hm ^^T^^i^W ((m, ti), xi; (ns, h),X2) 
T— >oo Uonj]^ 

5"-! r r ^wl/z+T2wl/Sh-s2W2/Sv 

dw2 I dwi— T- 



(5.20) 



(27ri)2 J- J- pwl/Z+nwl/Sh-siw^/S^ _ W2){Wi + W2) ' 



Consider now the a-dependent term, i^^^^ defined in fl3.20p . After the 
change of variable v = v — 1 and z = z — 1, (13.201) becomes 

d£.* dS* d.- (2i-l)(2i)-a) 



(27ri)3 /p^ _ /p^_^_, ^(^^ - i){z + v - l){w - v){w + v - a) 

The leading term again comes from the T~ ^/^-neighborhood of 1/2. After 
the change of variables 

w = I + |wiT-i/V5„ 5 = 1 + ^W2T-^^yS,, y = l + luT-^'^/S, (5.22) 

and controlhng the error terms as usual, we get 



lim ^^Ti/3i^(2)^(^^^^^)^^^.(^^^^^)^^^) ^ / /" dw;2 / dwi 

T^oo Conji (27ri)'^ Jr, Ai 



(5 23) 

(W2 - m)(w2 + m)(wi - m) V W - 'S'^'t y g«>»/3+Ti«;2/5,,-siii;i/S, ■ 

Finally, concerning the integration paths, from the local structure around the 
critical point, see Figure HI we obtain the conditions illustrated in Figure [5l 

□ 

There is still one region where the a and M dependence occurs in Figure[2l 
This is the transition between the Airy2 process and DBM. This is present 
for a e (1/2, 1) when n ~ (1 — a)H, or in terms of (6^, vr(6')), it occurs for 

2 — 2a + a^„ 
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Consider the scaling fl5.ll) with the condition fl5.24p and define the rescaled 
process as 

^ .Jt) - (t - 2yt(^3M))^ 

— X' ' 

In the large-T hmit Xt converges to the following limit process. 

Definition 21. The process ^dbm^2 is the process with m-point distributions 
at Ti < T2 < ■ ■ ■ < Tm given by the Fredholm determinant, 



P ( n {^DBM^2(rfc) < Sfc} ) = det (1 - XsKA^BM^,Xs)L2i{r,,...,r^}xR) 

(5.26) 



where Xs{'rk,x) = l{x > Sk). The kernel is defined by 

((s2-r|)-(.i-ri^))^ 



( \ exp ^ ^ 

^^DBM-,2l^l''5i; 7-2,52) = , = l[r2>ri] 

A/47r(r2 - Ti) 

Here 72 : e'^'/^oo — > e~^^^^oo, 71 : e~'^'^^f^oo e^'^^^^oo. Moreover, 71,72 pass 
on the left of and they do not cross. 

With this definition, let us state the result. 

Theorem 22. The process Xt defined in ( (5. 25\) converges to the process 
^dbm^2j niore precisely 

lim Xrir) = 5,^DBM^2(r/S;J (5.28) 

T — >oo 

in the sense of finite- dimensional distributions. The scaling coefficients Sh 
and are given by 

^^=((2-«)Tl-«)) ' ^^=l + vr'('^) + i^^'- ^'-''^ 

Proof of Theorem The first part of the proof is in complete analogy to 
the one of Theorem (THl with the main difference being that the critical point 
is at a instead of 1/2 (this explain why instead of n we get 0). We get the 
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following expression (with Si 

.(£2-£l)!' 



exp 



4(r2-Ti) 



+ 



1 



T2>T1 



du 



Using the identities 



1 

dW2 



+ 



72 



dwi- 



-- Si-Tt 
1 

(27ri)2 



72 



dwi 



71 



g-!«|/3+T2iy^-S2'!«2 ]_ 



UI2/3+T21d| — S2«J2 
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,«if/3+Ti«)f-si«;i - ■u)(Wl - U) 



Wi 



U 



M 



1 



1 



and 



{W2 — u){Wi — u) Wi — W2 \W2 — U Wi — U 
w ' ' ' ' 



\w\ 



w — u 



(5.30) 



(5.31) 



(5.32) 



n.>0 



we compute the integral over u, which has a simple pole when n = M — 1, 
letting to (15:271) . □ 

A priori one might want to modulate the slow particle rate like in (15. 2p 
but around some a instead of 1/2. This is however not a relevant change, 
since in a neighborhood of the curve in Figure [21 any point can be reached by 
fixing a and then choosing r to get the desired value of n/t or by choosing r 
and then modulating a. 

Up to a change in the time direction, the kernel K_^^^^,^^^ appeared in the 
context of sample covariance matrices [2] (for Ti = T2 = r), and the extended 
version in TASEP with step initial conditions [23], directed percolation with 
two set of parameters [10], Brownian Motions with outliers [1]. 

By looking at the diagram of Figure [2] it is quite apparent that one should 
have the following limits: 



lim A2^i.mA^ 
lim A2^i,mA'^ + ^ 

V—tOO 

lim lim ^2^i,m,k(^ ^ ^ 
lim A-i,Af,K(T - n 

re— >oo 

lim ^DBM-»2(l" + V 



^2^i(r), 
2i/3^i(r/22/3), 

A2{r), 

"^DBM^2(t), 

^2(r). 



(5.33) 
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6 Regions where the slow particles do not 
matter 

For completeness, we describe what happens in the region where the presence 
of slow particle is irrelevant. In the region where the density of particles is 
constant, the fluctuation of particles' positions are described asymptotically 
by the Airyi process. If the density of particles is decreasing (linearly in our 
case), then one has the Airy2 process, and in the transition region where 
the density changes from constant to linearly decreasing, the process is the 
Airy2^i process. The computations are essentially the same as in [8], but 
easily extended to the setting of space-like paths. The only difference is that 
one has to control the new term coming from K'^'^\ 

Introduce the scaling on space-hke paths described by a function 7r(^) 
with Itt'I < 1: 

t(r, T) = {Trie - tT-'/') + 9 - tT-'I') T, 
n(T, T) = M + [k{9 - tT-'/^) -{9- tT-'^^)] T. ^ ' ^ 

Case 1, n > max{^,l/4}t, i.e. 7r{9) > max{f^,|}e: The rescaled 
process 

. ^.(^)-(¥-2(n-M)) ^^ ^^ 

converges in the T ^ oo limit to the Airyi process, Ai, 

lim Xt(t) = S,Ai{r/Sh), (6.3) 

T — ''CO 

where Sy and are coefficients given by 

S, = {7:{9) + 9y/\ Sh = ^^^^SI (6.4) 

Case 2,ae (1/2, 1] and n G ((1 - af, l/4)t, i.e. ^9 > 7r{9) > ^^^9: 
The rescaled process 

y ( . Xn{t)-{t-2^t{n-M)) 

^T[r) = (6.5) 

converges in the large-T hmit to the Airy2 process, A2, 

lim Xt{t) = S,A2{r/Sh), (6.6) 

T— *oo 
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where Sy and Sh are coefficients given by 

-1/6 



2/3 



2(i-,f^Y' (6-7) 



Case 3, a E (1/2, 1] and n ~ t/4, i.e. n{9) = ^9: The rescaled process is 
given by 



' - (it-2(n-M)) 



Xr(r) = ( 



-Ti/3 



if n > t/4, 



Xt converges in the large-T hmit to the Airy2^i process, A2^i, 

hm Xt{t) = S,A2^i{T/Sh), (6.9) 
where Sy and Sh are coefficients given by 

7 Blocking wall regime 

In this section we study the case when M = oo and a = 2, so that the mean 
speed of the particles starting form 2N (called a-particles) is 1, which is equal 
to the jump rate of particles starting from 2Z_ (called normal particles). We 
want to describe the large time behavior of a finite number of normal par- 
ticles. For large time t, the a-particles fluctuate on a t^^^ scale, i.e., on a 
t^/^ scale their behavior is essentially deterministic, t^^^ is however the typi- 
cal scale of fluctuations of the normal particles, which perform random walks 
except for being blocked by their right-neighbor. Thus, one expects that for 
large time, our system should be related to a set of non-intersecting Brown- 
ian motions with some particular condition at the origin (like absorption or 
reflection). 

Theorem [2] stated in Section [2] is a direct consequence of the determi- 
nantal structure together with the following convergence of Koo (defined in 
Corollary [H]) to the kernel i^^GUE_ 
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Proposition 23. Let 

U = nt, Xi = U-U2tiY'\ (7.1) 

Then, 

t— >oo /i2 

w/iere = ln(7:j), with the conjugation factor = e~*'(trj/2)"*/^(— 2)"''t/^^, 
and ^^^g,^ j^j^ 



Before proving tlie result, let us present an integral representation of tlie 
antisymmetric GUE minor kernel, since it is in that form that we obtain the 
result. 

Lemma 24. The antisymmetric GUE minor kernel has the following integral 
representation (after conjugation) . Let Ti := e^' and e > 0. Then 



27ri ^[KA)^(n2,e.)] 



(27ri)2 /p^ e"'2^2-26V^«'2 - W2 + W2 / W2' 

(7.3) 

with the paths non-crossing, i.e., \w2\ < s, and Bi = 2"^e^'^"'"'"^^/^. 

Proof of Lemma\2^ We use the following two integral representations for the 
Hermite polynomials if„(x), 

Hn{x) = / dwe"" "^^'^u;", 

'V „ -^^^^ (7.4) 

as well as the identities (with < g < 1) which can be found in [24,27] 
1 ( {x-qyf\ ^Hk{x)Hk{y)q^ 

k=0 



/■oo 

/ dye-y"H4y)=e-^"Hn-iix), 

J X 

Hn{x) = {-lyHr^i-x). 



(7.5) 
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Then, for [ui, 6i) {n2, 62), we get (extending the sum to 00 because the 
extra terms are identically zero) that K''^^'^{{ni,ei),^i; (n2,^2),6) 

IS given 

by 



2"2 (27ri)2 z"^-- V^^e" 

Now, by the change of variables z = 1026^^^"^ = W2y/T2 and w = W\e^^^'^ 
W\Jt{ we obtain 



Performing the sum over 



'5 



/ \ 2€ 2 

E/ W2 \ UI2 r I I I I /^T c 

— = 2 < {7.t 

^ wi J wr — Wn 



£>1 



and replacing in (17. 7p one obtains fl7.3p . 

Now consider {ni,9i) -< (^12,62). Assume the following identity (proven 
below) 

oni^("l+l)/2 9 r— I- 1 
^ ^1 ^V'l / J„„„(Tl-r2W2_2(A,/?^-^2^/^)^« 



'2 

2 2 'H^^ g-{e2-eiK/2 



-T^e-^' E 2-^+1-^(^2 + ^^-^^ 

CO 



Then, the first two terms of (17.31) are equal to 

£=—00 

(7.10) 

and using the symmetry/antisymmetry properties of the Hermite polynomi- 
als, see (17. 5p . we get a zero contribution for all odd i. From this follows 
the result. It remains to show (17. 9p . We prove it by iteration, starting with 
rii = n2. The Gaussian integral gives 



9 '• 

20T 



g (n - T2 ) - 2 1 yrr- 6 



l-T2/n 



27ri JiK+g ■^7r(l — T-2/ri) 

e-g? ".^^ gn2+i-Kei)gn2+i-^(6)(r2/ri)"^+^-^/^ 
^/^ 2"2+i-^(n2 + l-^^' ^ ' 



-00 
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WW w 




Figure 7: Steep descent paths used for the different terms, with 
p= 1 — £:(t/2)~^/^. (a) is for term (1); (b) for term (2) and (3a); (c) for 
term (3b) in case Ti < T2; (d) for term (3b) in case Ti > T2. 



where we used (17. 5p with q = a/t2 / ri and replaced k hy n2 + 1 — i- Next, 
notice that the function F in (17.91) satisfies 

dyF^M = Fn,-i{x). (7.12) 

Thus, to get F„2 we have to integrate n2 — rii times F„2 . This is easily made 
using the integration formula in (17. 5p applied to e'^'^Hn^+i-ei^i) just n2 —rii 
times. This leads just to the shift in the index of the Hermite polynomials 
in dm]). □ 

Proof of Proposition We prove that under the scaling (17. ip 

lim ^(2ti)V2i^„o((ni,ti),a;i;(n2,t2),a;2) = (HD (7.13) 

t— >oo (^2 

where 0{t) = ln(r), with the conjugation factor Cj = e~*^(t/2)"*/^(— l)"'. 
The kernel Koo is given in Corollary [HI 

(1) Term coming from 0. We have 

-0«"-*^)-("-*^))(xi,a;2) = -— / —( ^fi^V'^'e^^'oM+t^/VW^ 

27ri/ro w \ w J 

(7.14) 

with 

^o(w) = /o,i(w) - /o,2(w), foA{w) = Ti{w -\n{w)), 

(7-15) 

gi{w) = /i,i(w) - /i,2(w), /i,i(w) = ^iV2Ti\n{w). 

Consider ti > r2 (the case ti = T2 and ni < n2 is pretty easy). As steep 
descent path we can use Vq = {w = pe"^, G (— 7r,7r]} for any p > 0, since 
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Re(w — ln(w)) = Re(w) — ln(p). The critical point of Qq is a.t w = 1 and there 
the Taylor series are 

go{w) = go{l) + ^in - T,){w - if + OHw - if), 



g,{w) = g^{l) + (Gv^ - 6v^) V2(^ - 1) + 0{{w - 1 



\2\ 



The leading contribution comes from the t~ ^/^-neighborhood of w = 1, which 
is however also a pole when n2> ni. Therefore we have to remain on its left 
and we choose p = I — (2/t)^/^e, e > 0. After controlling the error terms, we 
make the change of variable w = 1 + t^^^'^\/2W and take W = \K — e for any 
given £ > 0. So, the leading contribution of fl7.14p is given by 

V 1 27ri ■ ^ ' ' 



Finally, changing the variable W = —w and multiplication by Ci\/2ti/C2 
leads to the first term in (17.31) . 

(2) Term coming from the second term in ^3. 29\) . After the change of variable 
V = V ~ I we have 



(27ri)2 Jto i_- (M)"' e*/0.2(i')+ti/2/i,2(i)) yj(^v + yj- l){w - v) 

(7.18) 

The steep descent path for w is chosen as above, while for v we just take a 
circle around 1 of radius smaller than {2/tY^'^e. With the variables called 
v = l + t-^/'^V2V and w = 1 + t-^/^V2W, the leading contribution of (EH]) 
is 

' ' ■ ' ^ dV dW- 



C, V t (27ri)2 jTf^K, J,^_^ e.v-+2i,^,v 

(7.19) 

Finally, the change of variable W = —Wi, V = —W2 and multiplying by 
Ci\/2t7/C'2 leads to third term in (17. 3p (the part with 1/ (u'l — 102))- 

(3) Terms coming from the third term in Ii3. 29\) . Recall that we have to set 
a = 2 and after the change of variable v = v — I we obtain 



e 



dw d) dv- 



(27ri)2/r^ Jr,^_^ e^^^^Y^ w{w + v - 2){l + v - w)' 

(7.20) 
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(3a) Term coming from the pole at v = 1 of 1^7.20^ . In this case the sit- 
uation is almost identical as in case (2). The only difference is in the last 
factor, in particular, after the change of variables v = 1 + t''^^'^\/2V and 
w = l + t-^l^^piW the last factor goes to IjiyV + V) instead of IjiyV - V). 
After the final change of variable y&W = —Wi, V = —W2 and multiplication 
by Ciy^2ti/C2 leads to third term in (17.31) (the part with l/{wi + W2)). 

(3a ) Term coming from the pole at v = 2 — w of 1^7.20^ . This term reads 



(7.21) 
(7.22) 



where 

ho{w) = (ri + T2)w + r2(ln(2 - w) - 2) - n ln{w) 
hi (w) = V2n^i \n{w) - V2t^6 ln(2 - w) . 
There are two critical point of ho, namely 



= 1, uj2= , both in [0, 1]. (7.23) 

n + T2 

The steep descent path passes by the critical point the closest to the origin. 
For Ti < T2, we have lj2 < ^ and the steep descent analysis gives readily a 
contribution of order 

C 

^^tho{i^2)-thQ(i) ^ ^7_24) 
Ci 

It is easy to see that, with := T2/T1, 

/io(^2)-/io(l) = ri(l-;u)(l + ln(2)-ln(/i + l)) < 0, for all /i > 1. (7.25) 

Thus in the t —>■ 00 limit, the contribution goes to zero exponentially fast for 
Ti < T2. Finally, consider ti > T2. Then, 1 < uj2- We choose the path as in 
case (1), but this time the Taylor series give 

ho{w) = ho{l)+^^^^^{w~ir, ^^26) 
hi{w) = /io(l) + i^iVn+ ^2V^)V2{w - 1) + 0{{w - 1)2). 

Also, we have a different sign in the pref actor and a factor (—1)"^ in the 
term {w — l)"i~^"2^ This leads to the second term in (17. 3p and explains the 
differences with the first term of (17. 3p . namely the (— l)"2+i the change 

6 - -6- □ 
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